Topological phases of Bi(lll) bilayer in an external exchange field 
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Using first principles methods, we investigate topological phase transitions as a function of ex- 
change field in a Bi(lll) bilayer. Evaluation of the spin Chern number for different magnitudes 
of the exchange field reveals that when the time reversal symmetry is broken by a small exchange 
field, the system enters the time-reversal broken topological insulator phase, introduced by Yang 
et al. in Phys. Rev. Lett. 107, 066602 (2011). After a metallic phase in the intermediate region, 
the quantum anomalous Hall phase with non-zero Chern number emerges at a sufficiently large 
exchange field. We analyze the phase diagram from the viewpoint of the evolution of the electronic 
structure, edge states and transport properties, and demonstrate that different topological phases 
can be distinguished by the spin-polarization of the edge states as well as spin or charge transverse 
conductivity. 
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I. INTRODUCTION 

Topological insulators (TIs) have drawn quite in- 
tensive attention recently owing to the gapless sur- 
face/edge states, topologically protected against pertur- 
bations which do not close the bulk gap>ii^ Similarly 
to graphene^ exciting physical properties are expected 
if the TI surface/edge states are made gapped, result- 
ing in e.g. half-quantized surface Hall effect^ topological 
magneto-electric effectSf^i^ and Majorana zero states<i It 
is suggested that such effects can be triggered via the 
proximity effect of TIs to magnetically ordered materials 
or superconductors, in which case new terms which im- 
pose or break certain symmetries are introduced into the 
effective Hamiltonian. 

One particularly interesting situation arises if the time- 
reversal symmetry (T) in a TI is broken via a controlled 
perturbation. The definition of the Z2 topological in- 
variant, used to classify two-dimensional (2D) insulators 
into ordinary insulators and TIs, is hinged on the pres- 
ence of the time- reversal symmetry in the system<^ It still 
remains an open question of how to generalize the Z2 
topological invariant when 7" is broken. In such a case, it 
is suggested^ that it is plausible to use the so-called spin 
Chern numberii^ 



C. = \{C, 



C-), 



(1) 



where C+ and C_ are the Chern numbers for the spin " up" 
and "down" manifolds of the occupied states. Since the 
spin-orbit coupling (SOC) in general induces spin mix- 
ing, the spin Chern number is well defined only when the 
system exhibits a gap in the spectrum of ct^ in addition 
to being a band insulator. In this case Cg proves to be 
equivalent to the Z2 number when time-reversal symme- 
try is present iii For models of TIs^^ and graphenc;^ the 
spin Chern number picture proved to be valid, while a 
more detailed analysis of this concept in real materials 
investigated from first principles is still lacking. 

Owing to T symmetry, the surface/edge states of TIs 
are spin degenerate, making it hard to manipulate and 



exploit their transport properties in experiments based 
on magnetic detection. While the quantum Hall effect 
usually requires a very large magnetic field to emerge, 
in this sense, the quantum anomalous Hall (QAH) effect 
in fcrromagnets would be much more suited for future 
spintronics. The QAH effect, which has attracted a lot 
of attention recently but so far has not been verified ex- 
perimentally, occurs in 2D Chern insulatorsi^ with bro- 
ken 7~-symmetry, exhibiting overall spin-polarized edge 
states carrying quantized electric charge. Up to now, it 
is believed that TIs can be turned into Chern insulators 
via magnetic dopingi^J^ but it still remains an exper- 
imentally unsolved problem despite recent advances in 
this direction ji^ In this light, it is necessary to explore 
how tunable 2D TIs are under T-brokcn perturbations, 
in particular, exchange interactions. 

Bismuth is the heaviest atom with a stable isotope 
and strong spin-orbit coupling makes it an important in- 
gredient of the newly discovered TIs such as Bi2Se3i^ 
and Bii_^Sba;.— Ultrathin Bi(lll) films, which can be 
produced experimentally on different substrates , ^^i^° are 
also predicted to be 2D TIsi^' 



1^ In our study we take 
Bi(lll) bilayer as a representative of the latter class of 
TIs. In this work, we investigate the topological phase 
transitions in a 2D TI - the Bi(lll) bilayer— — in an ex- 
ternal exchange field using first principles methods. We 
found that at small exchange fields the bilayer is in the 
time-reversal broken topological insulator phase, which 
can be characterized by non-trivial values of spin Chern 
numbers, transverse spin Hall conductivity, and spin po- 
larization of the edge states. Increasing the magnitude 
of the exchange field further drives the system first into 
a metal and then into a QAH phase. Detailed analy- 
sis of the electronic structures reveals that spin-mixing 
is important for closing the bulk gap and Chern num- 
ber exchange between valence and conduction bands. At 
last, we analyze the spin-polarization of the metallic edge 
states in different phases, show that its sign can be con- 
trolled by the strength of the exchange field, and prove 
that its magnitude is large enough to be observed in scan- 
ning tunneling microscopy experiments. 



II. METHOD 

In a Bi(lll) bilaycr two layers of Bi atoms form a hon- 
eycomb lattice when projected onto the plane of the film. 
The relaxed bulk in-plane lattice constant and the dis- 
tance between the two layers were 4.52 A and 1.67 A, re- 
spectively (see also Fig. 4). Our theoretical investigations 
are based on density functional theory.— We apply the lo- 
cal density approximation^^ to the exchange-correlation 
potential and use the full-potential linearized augmented 
plane wave method (FLAPW) as implemented in the 
FLEUR codci^ The self-consistent calculations with SOC 
were carried out with the cut-off parameter k,nax of 
3.8 bohr~^ and 50 fc-points in the full two-dimensional 
Brillouin zone (BZ). The muffin-tin radius of 2.5 bohr 
was used. The band structure of the system with SOC 
is presented in Fig. 2(f). The Wannier functions tech- 
nique was used on top of self-consistent DFT calcula- 
tions to derive an accurate tight-binding Hamiltonian 
of the systemi^"— Moreover, we have constructed the 
matrix elements {ipmk \ <^a \ V'nk) of the Pauli matrices 
(Ja{a — X, y, z), where V'nk are the occupied Bloch wave 
functions. These matrices were evaluated from the DFT 
calculations and transformed into the real space repre- 
sentation in terms of maximally localized Wannier func- 
tions (MLWFs)i^ This allowed us to take into account 
an exchange field applied perpendicularly to the surface 
of the bilayer by adding a. a^ ■ B term on top of the 
original Hamiltonian, and to provide an ability to ac- 
curately evaluate spin polarization Pa (see caption of 
Fig. 4 for details). In general, applying an exchange 
field lifts the spin degeneracy, where spin-up (spin-down) 
states get higher (lower) in energy spectra, and the differ- 
ence of the occupation of spin-up and spin down leads to 
the spin-polarization Pa which can be positive (negative) 
along the direction of exchange fields. In Fig. 2(a)-(e), 
Fig. 2(k)--(o), and Fig. 4 shown below, spin-up (minor- 
ity) and spin-down (majority) states are indicated in red 
and blue, respectively. 

For a band insulator the transverse anomalous Hall 
conductivity (AHC) is ^ times the conventional integer 
(first) Chern number, where^i 

-. OCC r, 

c = c+ + c- = -^Y.J d'kniy{-k), (2) 
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FIG. 1: Phase diagram of Bi(lll) bilayer with respect to the 
exchange field B. Blue squares (red circles) mark the band 
gap (anomalous Hall conductivity AHC) as a function of B. 
Regions of different topological phases, namely, T-conserved 
TI (TRS TI), T-broken TI (TRB TI), metal, and Chern insu- 
lator (QAH) are shaded difi'erently and labeled accordingly. 



evaluate the spin Hall conductivity, we have replaced v^ 
with the spin-velocity operator sf, = {vu,crz} in Eq. (3). 
To calculate the spin Chern number, we followed the 
procedure of Ref. [l3|. First, the CTz matrix for the occu- 
pied states {iprn'k \ Cz | V'nk) {"m, n are the band indices of 
occupied states at each k) is constructed and diagonal- 
ized. Note that due to the spin- flip part of the spin-orbit 
coupling^, the eigenvalues of Cz matrix are not necessar- 
ily ±1. Then, using the "spin- up" or "spin-down" eigen- 
vectors \(j)^) of the (Jz matrix, the occupied states are 
projected into the "spin- up" and "spin-down" manifolds 
■0^ = (0^(k) I '0(k)), where |'0(k)) denotes the eigen- 
vector of the occupied states. At last, the spin Chern 
number is evaluated by integration of the Berry curva- 
ture for each manifold where the derivatives of the wave 
functions ip with respect to k are obtained using the 
finite difference methodsi^ For the spin Chern number, 
it is only well defined if the spectrum of CTz is gappedJ^ 
In our calculations, we observed that in the range of ex- 
change field considered (0 eV < B < I eV) the spectrum 
of (Tz remains finite which justifies the usage of spin Chern 
numbers in this work. 



in which the summation goes over all occupied states and 
0„(k) is the Berry curvature of the n-th band, given by 
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with w„k as the lattice-periodic part of ijjnk with energy 
eigenvalues e„k, i'v {v = x, y) is the velocity operator. To 



III. CHERN NUMBERS AND THE PHASE 
DIAGRAM 

Here, we calculate the Chern numbers and analyze the 
topological phases of Bi(lll) bilaycr as a function of the 
externally added exchange field B, presenting the phase 
diagram, calculated band gap and AHC in Fig. I. When 
5 = 0, the system is a 2D TI with Z2 number of 1 (non- 
trivial). The Chern numbers for the spin- up (minority) 
and spin-down (majority) manifolds are C± = ^1, re- 
spectively, leading to a spin Chern number Cs — —1 
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FIG. 2: Electronic structure of Bi(lll) bilayer in an external exchange field perpendicular to the surface, (a)-(e) display the 
2D "bulk" band structure for an exchange field of 0.0, 0.15, 0.35, 0.5, and 0.8 eV (indicated at the top), respectively. Red 
(blue) color stands for positive (negative) spin polarization along the field of the states in arbitrary units. The band structures 
are plotted along the K-F-K' path in the vicinity of the F point. In (a) the spin- up states are shifted by 10 meV upwards in 
energy to make the degeneracy of the states more visible. In (f) the 2D "bulk" band structure of Bi(lll) bilayer is displayed. 
Each band is doubly degenerate, (g)-(j) display the distribution of the Berry curvature of occupied states in the BZ (to a minus 
sign) with exchange field marked on the top. Blue (red) color indicates positive (negative) values of the Berry curvature (in 
a.u.^). The evolution of the band structure of an 80-atom wide zigzag ribbon is shown in (k)-(o). Edge states are colored with 
the expectation value of the a^, rendered in blue and red for negative or positive values, respectively. Size of points indicates 
the weight from atoms located on the upper edge of the ribbon. The gray shaded regions denote the projected 2D "bulk" band 
structure from 2D BZ onto the ID k vector of the ribbon. In (k) the spin-up states are shifted by 10 meV upwards in energy 
to make the degeneracy of the states more visible. Horizontal dotted lines in (a)-(e) and (k)-(o) indicate the constant energies 
(chosen to be in the bulk gap in insulating phases) for which the spin polarization is analyzed as shown in Fig. 4. 



(cf. Eq. (1)) in agreement with previous ealculationsj^ In 
this case, the spin Chern number is equivalent to the Z2 
number. Introducing and increasing the exchange field 
breaks the T'-symmetry and causes an exchange splitting 
between the spin-up and spin-down valence and conduc- 
tion bands, bringing thus the minority valence bands to- 
wards majority conduction bands,— Fig. 2(b). Further 
increase of the exchange field leads to a closure of the gap 
aX B = 0.31 cV, Fig. 2(c). As evident from Fig. 2(a)-(c), 
until the bulk gap is closed, the original inverted band 
structure at the F-point remains topologically nontrivial. 
This is confirmed by calculating the spin Chern number 
C± = =Fl for < B < 0.31 eV, identifying the topologi- 
cal phase of Bi(lll) bilayer in this range of the exchange 
field as the T-broken, or, TRB TI phase (cf. Fig. 1), 
analogous to that considered in graphcne^ 

To verify that the T-brokcn phase could display non- 



trivial topological transport properties, we calculate the 
transverse spin Hall conductivity (SHC) as a function of 
the spin-orbit strength in the system. A, and present the 
results in Fig. 3. Consider first the bilayer at B = 0. 
At small SOC strength the system is in a trivial insula- 
tor phase, which can be characterized with a zero SHC. 
Upon increasing A, the bilayer first goes through a metal- 
lic phase, and at around A = 0.7 of the SOC strength of 
Bi atoms a transition to the TI phase with SOC-driven 
band inversion occurs, which is accompanied by a non- 
zero SHC of the magnitude of around —0.7^3^ The de- 
viation from the quantized value of — ^ is due to the non- 
vanishing spin-non-conserving part of the spin-orbit in- 
teraction: artificially switching it off in our calculations^ 
we indeed acquire a quantized value of — ^ for the SHC 
(dashed line in Fig. 3), in accordance to the original sce- 
nario proposed first for HgTc/(Hg,Cd)Te quantum spin 
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FIG. 3: Spin Hall conductivity (SHC) and anomalous Hall 
conductivity (AHC) of Bi(lll) bilayer with respect to the 
scaled strength of the spin-orbit coupling of Bi atoms A. A = 1 
corresponds to the unsealed SOC strength. Black circles (red 
squares) denote the SHC with B = 0.0 eV {B = 0.15 eV), 
while purple diamonds mark the AHC with B — 0.15 eV. 
Horizontal dashed line at SHC=-^ stand for the quantized 
magnitude of SHC if the spin-flip part of SOC is switched off. 



Hall hctcrostructureS)^"— . By looking at the SHC for 
the Bi bilayer at B = 0.15 eV, we observe an almost 
identical behavior to that vifithout the field, with only 
modifieations in the boundaries of the metallie region, 
SHC in the metallic regime, and a slight difference in 
the value of the SHC at large A. Again, by neglecting 
the spin-non-conserving part of SOC, we restore a quan- 
tized value of -^ for the SHC in the TRB TI phase, 
which clearly manifests a topologically non-trivial phase, 
as far as the topological transport is concerned. Note 
that for both phases, TI and TRB TI, the transverse 
charge conductivity in exactly zero when A lies outside 
of the metallic region. The ability to characterize the T- 
brokcn TI phase with observable in an experiment non- 
zero spin Hall conductivity could be particularly impor- 
tant for distinguishing this phase from a trivial insulator, 
since, as we shall see in the next section, the topological 
protection of the edge states in the system for i? 7^ is 
lifted. 

Returning now to the phase diagram in Fig. 1, we ob- 
serve that the size of the band gap decreases linearly as 
the exchange field is increased, and for B in between 0.31 
and 0.42 eV the bulk gap is closed with Bi(lll) bilayer 
in the metallic phase, see Fig. 1. However, increasing the 
exchange field further reopens the gap, which results in 
the occurrence of a topologically non-trivial phase. For 
B > 0.42 eV the spin Chern number of the occupied ma- 
jority states reverses sign, leading to a zero spin Chern 
number Cs , while the Chern number C acquires a value of 
—2. That is, in this range of the exchange field the sys- 
tem is in the QAH, or, Chern insulator phase. It is worth 
to recall here that in case of a metal in which the spin- 
non-conserving part of the spin-orbit interaction is hy- 



pothctically switched off^^ the intrinsic anomalous Hall 
conductivity is given by the sum of the Hall conductivi- 
ties for spin- up and spin-down bands, while the spin Hall 
conductivity is given by their difference (spin is a good 
quantum number in such a case)i^2i^ Evidently, such a 
correspondence between the Chern number and the spin 
Chern number, i.e., their representation as the sum and 
the difference of the Chern numbers for spin up and spin 
down bands, holds true for Bi(lll) bilayer with gapped 
spectrum of Uz, based on the argument that a unitary 
transformation (eigenvectors of az in our case) of the oc- 
cupied states will keep the Chern number invariant 1^ 

Examining the electronic structure of the system, pre- 
sented in Fig. 2, reveals that it is the spin- mixing accom- 
panied by the exchange of the Chern number between 
valence and conduction bands that leads to the phase 
transition from the T-conserved TI phase to the QAH 
phase through the T-broken TI and metallic phases. The 
occurrence of the metallic phase has not been predicted 
previously e.g. for grapheno^i or Bi2Se3ji^ It emerges due 
to an overlap in energy of bands with opposite spin char- 
acter in the vicinity of the F-point (cf. Fig. 2(c)), which 
is characteristic of a material with a non- monotonous dis- 
persion of the bands on both sides of the F-point. Such 
a metallic phase cannot be achieved e.g. for pure Dirac 
bands, which have a local maximum (minimum) at the 
high symmetry (in this case) F-point, as is evident from 
Fig. 2(c). 

Our transport calculations according to Eq. (3) show 
that a finite AHC develops in the metallic and QAH 
phases, see Fig. 2. Microscopically, the non- vanishing 
AHC can be traced back to the development of the sin- 
gular Berry curvature in k space as the _B-field is in- 
creased, see Fig. 2(g)-(j). For B below 0.31 eV the Berry 
curvature has a non-trivial distribution in the BZ with 
the regions of positive and negative values, comparatively 
small in magnitude (Fig. 2(g)), while the BZ integral of 
the Berry curvature amounts to zero. 

At the onset of the metallic phase there is a noticeable 
hybridization between the conduction and valence bands 
of opposite spin, which gives rise to a very large negative 
contribution to the Berry curvature around the points in 
the BZ at the Fermi level (Ep) where the orbital charac- 
ter and the spin character of the valence and conduction 
bands is changing, or, in terminology of Ref. [331 j along 
a "hot loop" in the BZ, Fig. 2(c)-(e) and (h)-(j). In such 
a situation, the Berry curvature and the AHC, propor- 
tional to the integrated value of the Berry curvature over 
the BZ, are dominated by the spin-flip transitions me- 
diated by the non-spin-conserving part of the spin-orbit 
interactioni^ At the F-point, the sign of the Berry cur- 
vature reverses for B > 0.3 eV, accompanying the switch 
in parity of valence and conduction bands. 

In the metallic phase, the AHC as a function of B is 
monotonously increasing (Fig. 1), while the radius of the 
hot loop is increasing owing to the shift of the hybridiza- 
tion point between the valence and conductance bands 
away from F (Fig. 2). When the gap is reopened for 



B > 0.42 eV, the AHC reaches a value of -2^, with the 
Chern number C = —2, see Fig. 2. In the QAH phase, 
we observe from our calculations that out of six valence 
bands, the topmost two bands dominate the contribution 
to the AHC for exchange fields up to 0.60 cV. When the 
exchange field is increased further, deeper spin-up bands 
will be pushed closer to the Fermi energy and take the 
spin-down character derived from the conduction bands 
(cf. Fig. 1(c)), also providing contribution to the Berry 
curvature. 



IV. EDGE STATES 

We turn now to the properties of the edge states in 
Bi(lll) bilaycr zig-zag ribbon with the width of 80 atoms 
in y direction and periodic (infinite) in x dircctioui^ The 
tight-binding Hamiltonian of the ribbon was constructed 
in terms of the MLWFs by a mere termination of the 2D 
(infinite) Hamiltonian after 40 unit cells in the y direc- 
tion. For Bi(lll) ribbon, when T symmetry is not bro- 
ken, there are twelve edge states in total, as evident from 
Fig. 2(k). Six of them are located on each side of the rib- 
bon, while out of these six there are three "spin- up" left- 
movers and three "spin-down" right- movers. The axis of 
spin polarization of the edge states is determined by the 
details of the electronic structure and spin-orbit interac- 
tion, and can be also energy-dependent 3^ 

Applying the exchange field to the ribbon breaks the 
7~-symmetry, lifts the spin degeneracy and polarizes the 
edge states along the direction of the applied field, see 
Fig. 2(l)-(o). A close inspection of the edge bands in 
the case of the TRB TI phase. Fig. 2(1) reveals that the 
edge states do not connect the valence and conduction 
bands in the system, as in the case of the edge states 
of a conventional TI, cf. Fig. 2(k). This means that 
the Hamiltonian of the system can be transformed such 
that the spectrum at the edges exhibits a gap at the 
Fermi energy. Moreover, owing to the breaking of the 
7~-symmetry, back-scattering is in general allowed which 
results in the dissipative transport properties of the edge 
states.— Noticeably, in the QAHE phase. Fig. 2(n,o) the 
connection of the valence and conduction bands of the 
bilayer via the edge states is restored again, resulting 
in a dissipationless longitudinal transport of the ribbon, 
realized by the topologically protected edge states. 

Of particular interest here is the spatial distribution 
of the spin-polarization of the edge states around the 
Fermi energy, P^, which plays a very important role in 
spin-polarized scanning-tunneling microscopy (SP-STM) 
experiments!^ We calculate this spin polarization as a 
function of the _B-ficld in real space across the ribbon, 
and present the results in Fig. 4. Let us consider first 
the case of a very small exchange field B = 0.05 eV and 
look at the spin-polarization at Ep across the ribbon 
(Fig. 4). We observe that the spin polarization is: (i) 
rather localized at the edges of the ribbon, (ii) reaches 
very large values, (iii) displays oscillations in sign across 
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FIG. 4: (color online.) Spatial distribution of the spin polar- 
ization Pz at the Fermi energy of the edge states in a zigzag 
Bi(lll) ribbon with tlie width of 80 atoms. P^ is obtained as 
a sum of the expectation values of a^ for all edge states at the 
constant energies (marked with dashed lines in Fig. 2(a)-(e)), 
multiplied by the weight of the corresponding wavefunction 
on given atom. The structure of the edge states on each side 
of the ribbon is marked with circles to the left (left edge) and 
to the right (right edge) of the value of the corresponding 
exchange field B. While the direction of propagation of the 
edge states is marked with white dots (white crosses) for the 
direction towards the reader (away from the reader), the red 
(blue) color of the circles marks predominantly spin-up (spin- 
down) character of the edge states. Lower inset depicts the 
distribution of P, for B = 0.05 eV (in arbitrary units with 
red for positive and blue for negative values) along and across 
the ribbon, where yellow spheres stand for the positions of Bi 
atoms. 



the ribbon. These properties can be understood by refer- 
ring to the electronic structure of the edge states for this 
value of the field, which is quite similar to that depicted 
in Fig. 2(1), only with the smaller separation between 
the upper (predominantly spin-up) and lower (predomi- 
nantly spin-down) subgroups of two edge states, so that 
there are twelve states at Ep all together. The direction 
and spin character of the edge states at Ep is sketched 
in Fig. 4. We observe that owing to the presence of the 
non-zero exchange field and breaking of T symmetry, the 
edge-states get spin-polarized, while on each side of the 
ribbon the number of right and left movers is equal for 
B < 0.42. This results in zero transverse charge conduc- 
tance, in accordance to the phase diagram. Fig. 1. 

For larger wavevectors \kx • a| « tt, see Fig. 2(1), the 
upper group of two states has a strong spin-up character, 
and for B = 0.05 eV (not shown) it crosses the Fermi 
energy, giving rise to a strong spin-polarization in the 



direction of the _B-field. Since the states with larger k^ 
(away from the bulk bands) are more localized at the 
edges of the ribbon than the states with smaller k^ (closer 
to the bulk bands) ^21 positive P^ in Fig. 4 for this field 
is localized at the ribbon edges. Correspondingly, the 
states at Ep close to F are the source of negative spin 
polarization, which is smeared more across the ribbon 
owing to stronger hybridization with the bulk states, thus 
giving rise to the change in the sign of Pz away from the 
ribbon edges. We speculate that such oscillatory spatial 
dependence of P^ should be a common phenomenon in a 
situation where the electronic structure of the edge states 
of a 2D TI exposed to an exchange field is non-trivial. 
The spatial distribution of the Pz in the unit cell for this 
case is illustrated in the inset of Fig. 4. 

Increasing the exchange field further pushes the tail of 
long wave length spin-up edge states out of the bulk gap 
completely around B = 0.15 eV, Fig. 2(1). In this case, 
there arc only 8 edge states left, whose direction of the 
spin-polarization at the Fermi energy is opposite to the 
direction of the i?- field. Since the Bi(lll) bilayer exhibits 
a T-brokcn TI phase for this value of B, there are two 
left and two right movers on each side of the ribbon, with 
zero transported charge, Fig. 4. In correspondence to the 
predominantly spin-down character of the edge states in 
k space, the Pz is large and negative in the vicinity of 
the edge of the ribbon. Noticeably, the spin polarization 
decays much slower with the distance towards the cen- 
ter of the ribbon, when compared to the case considered 
previously, owing to the proximity of the edge states in 
k space to the bulk states. We remark here that the pre- 
dicted change in sign of Pz is an interesting phenomenon, 
which can be exploited experimentally e.g. via tuning the 
strength of the exchange field by deposition of different 
adatomai^i^ (see also comments in section V). 

The evolution of the edge states for 0.15 < B < 0.80 
eV can be seen in Fig. 2(l)-(o). As the exchange field is 
increased, the magnitude of Pz gradually decreases, since 
the conduction and valence bands become very close in 
energy, and the kp vectors of the edge states approach 
the r point. This is particularly clear for the edge states 
in the QAH phase at B ~ 0.8 eV, for which Pz and its 
decay rate into the ribbon are very small. Fig. 4. For 
the regime of B in between 0.42 and 0.7 eV, i.e., at the 
onset of the QAH phase, the number of edge states is the 
same as that in the T-broken TI phase at B = 0.15 eV, 
but the emergent quantized charge conductance of the 
edge states proportional to C on each side of the ribbon 
is evident, Fig. 4. When i? > 0.7 cV four more edge 
states arc pushed out of the bulk gap window, and only 
four edge states crossing Ep are left overall. In this case, 
we have a situation of two right/left movers on each side 
of the ribbon, with a small positive spin polarization. 
Fig. 4. The microscopic mechanism for spin-polarized 
electron conduction in these edge states without possi- 
bility to back-scatter remains and open and interesting 
topic for future studies. 



V. SUMMARY 

In this work, based on first principles calculations, 
we have considered the emerging topological phases of 
Bi(lll) bilayer in an external exchange field B. We iden- 
tify four different phases as the exchange field is varied: 
TI phase for i? = 0, T-broken TI phase for small fields, 
metallic phase for intermediate fields, and Chern insula- 
tor, or, QAH phase for large exchange fields. Wc con- 
sistently identify each of the insulating phases in terms 
of the Chcrn number C, spin Chern number C^, spin and 
charge transverse conductivity. We attribute the origin of 
phase transitions to the spin-orbit mediated spin-mixing 
between valence and conduction bands whose topology 
in k space is controlled by the applied field. 

In particular, wc focus on the electronic structure and 
development of the edge states in a zigzag ribbon of 
Bi(lll) bilayer as a function of the S-field. We show 
that the spin polarization of the edge states along the 
direction of the exchange field at the Fermi level can be 
significant in magnitude, and its sign, as well as its spa- 
tial distribution can be controlled by changing the mag- 
nitude of the applied field. Since reversing the direction 
of the exchange field leads to reversed propagation direc- 
tion and reversed spin polarization of the edge states, the 
exchange field provides a tool to tune the properties of 
one-dimensional spin-polarized transport arising at the 
edges of insulators with non-trivial topological proper- 
ties. We also speculate that the spin polarization of 
the edge states, which characterizes different topologi- 
cal phases, can be observed in e.g. scanning tunneling 
microscopy experiments. 

We remark that the experimental realization of the 
tunable exchange field applied to a thin Bi(lll) bilayer 
can be sought in two directions. One way lies in finding 
a suitable insulating substrate for deposition of Bi(lll) 
which exhibits essential magnetism at least at the inter- 
face with the bilayer. In this case the direction of the B- 
field can be easily (depending on the magnetocrystalline 
anisotropy energy) reversed, while the magnitude of the 
exchange field can be tuned via e.g. the magnctoclcctic 
effect if the substrate used for deposition has multifcrroic 
properties. Second route lies in deposition of transition 
metal adatoms on the Bi(lll) surface. Recently, it has 
been demonstrated theoretically that such adatom de- 
position can induce considerable exchange splitting with 
magnitude as large as several tenth of eV on the topo- 
logical insulator underneath;^ and for Bi2Se3, the ex- 
change field is large enough to get the system into the 
QAH phascJ^. The magnitude of the exchange field in 
this case can be controlled by an appropriate choice of 
the transition metal, or by deposition of the system on 
an insulating non-magnetic substrate, which can cause 
essential modifications in the distance between the TI 
and deposited adatoms. We also note, that, generally 
speaking, the magnitude of the exchange field necessary 
to achieve the QAH phase in the phase diagram of Fig. 1 
depends directly on the initial value of the gap in the TI 



phase without the field. 



VI. ACKNOWLEDGMENTS 



This work was supported by the HGF-YIG Programme 
VH-NG-513 and by the DFG through Research Unit 912 
and grant HE3292/7-1. Computational resources were 
provided by the Jiilich Supercomputing Centre. 



We acknowledge very helpful discussions with Emil 
Prodan, Sheng Li, Klaus Koepernik and Marjana Lezaic. 



* corresp. author: h.zhang@fz-juelich.de 

^ M.Z. Hasan, and C.L. Kane, Rev. Mod. Phys. 82, 3045- 

3067 (2010). 
^ X.-L. Qi, and S.-C. Zhang, Rev. Mod. Phys. 83, 1057-1110 

(2011). 
^ A.H.Castro Neto, F. Guinea, N.M.R. Peres, K.S. 

Novoselov, and A.K. Geim, Rev. Mod. Phys. 81, 109-162 

(2009). 
" X.-L. Qi, T.L. Hughes, and S.-C. Zhang, Phys. Rev. B 78, 

195424 (2008). 
^ X.-L. Qi, R. Li, J. Zang, and S.-C. Zhang, Science 323, 

1184-1187 (2009). 
® A.M. Essin, J.E. Moore, and D. Vanderbih, Phys. Rev. 

Lett. 102, 146805 (2009). 
'' L. Fu, and C.L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 

* L. Fu, and C.L. Kane, Phys. Rev. B 76, 045302 (2007). 

^ Y. Yang, Z. Xu, L. Sheng, B. Wang, D.Y. Xing, and D.N. 

Sheng, Phys. Rev. Lett. 107, 066602 (2011). 
^° E. Prodan, Phys. Rev. B 83, 195119 (2011). 
" E. Prodan, New. J. Phys. 12, 065003 (2010). 
^2 C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, 

Phys. Rev. Lett. 101, 146802 (2008). 
" R. Yu, W. Zliang, H.-J. Zhang, S.-C. Zhang, X. Dai, and 

Z. Fang, Science 329, 61-64 (2010). 
" C.-Z. Chang, J.-S. Zhang, M.-H. Liu, Z.-C. Zhang, X. 

Feng, K. Li, L.-L. Wang, X. Chen, X. Dai, Z. Fang, X.-L. 

Qi, S.-C. Zhang, Y. Wang, K. He, X.-C. Ma, Q.-K. Xue, 

laxXiv: 1108.4754. 
^^ H. Jin, J. Im, and A.J. Freeman, Phys. Rev. B 84, 134408 

(2011). 
^^ F.D.M. Haldane, Phys. Rev. Lett. 61, 2015-2018 (1988). 
" K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. 

Katsnelson, LV. Grigorieva, GV. Dubonos, and A. A. 

Firsov, Nature 438, 197-200 (2005). 
^* Y. Xia, D. Qian, D. Hsieh, L. Wray A. Pal, H. Lin, A. 

Bansil, D. Grauer, Y.S. Hor, R.J. Cava, and M.Z. Hasan, 

Nature Physics 5, 398-402 (2009). 
^^ T. Nagao, J.T. Sadowski, M.Saito, et ai, Phys. Rev. Lett. 

93, 105501 (2004). 
^^ T. Hirahara, G. Bihlmayer, Y.Sakamoto, et ai, Phys. Rev. 

Lett. 107, 166801 (2011). 
^^ M. Wada, S. Murakami, F. Freimuth, and G. Bihlmayer, 

Phys. Rev. B 83, 121310 (R) (2011). 



22 S. Murakami, Phys. Rev. Lett. 97, 236805 (2006). 

2^ D. Xiao, W. Zhu, Y. Ran, N. Nagaosa, and S. Okamoto, 

arxiv:org:1106.4296. 
2* Z. Liu, C.-X. Liu, Y.-S. Wu, W.-H. Duan, F. Liu, and J. 

Wu, Phys. Rev. Lett. 107, 136805 (2011). 
2^ P. Hohenberg, and W. Kohn, Phys. Rev. 136, B864-B871 

(1964). 
2^ V.L. Moruzzi, J.F. Janak, A.R. Williams, Calculated Prop- 
erties of Metals, New York: Pergamon, (1978). 
2^ For the description of the code see www.flapw.de 
2* X. Wang, J.R. Yates, I. Souza, and D. Vanderbilt, Phys. 

Rev. B 74, 195118 (2006). 
2^ F. Freimuth, Y. Mokrousov, D. Wortmann, S. Heinze, and 

S. Bliigel, Phys. Rev. B 78, 035120 (2008). 
^° A.A. Mostofi, J.R. Yates, Y.-S. Lee, I. Souza, D. Vander- 
bilt, and N. Marzari, Comput. Phys. Commun. 178, 685 

(2008). 
^^ F.D.M. Haldane, Phys. Rev. Lett. 93, 206602 (2004). 
^2 This is in contrast to the case considered in Ref. [1^ , where 

states of the same spin character from the conduction and 

valence bands are pulled to the Fermi energy by increasing 

the spin-orbit strength. 
^^ L Souza, N. Marzari, and D. Vanderbih, Phys. Rev. B 65, 

035109 (2001). 
^"^ S. Murakami, N. Nagaosa, and S.-C. Zhang, Phys. Rev. 

Lett. 93, 156804 (2004). 
^^ B.A. Bernevig, T.L. Hughes, and S.-C. Zhang, Science 

314, 1757 (2006). 
^•^ M.Konig, S. Wiedmann, C. Briine, et al, Science 318, 766 

(2007). 
^■^ L. Sheng, D.N. Sheng, C.S. Ting, and F.D.M. Haldane, 

Phys. Rev. Lett. 95, 136602 (2005). 
^* H. Zhang, F. Freimuth, S. Bliigel, Y. Mokrousov, and I. 

Souza, Phys. Rev. Lett. 106, 117202 (2011). 
^^ R. Resta, Rev. Mod. Phys. 66, 899 (1994). 
*° C. Bruene, A. Roth, H. Buhmann, E.M. Hankiewicz, 

L.W. Molenkamp, J. Maciejko, X.-L. Qi, and S.-C. Zhang, 

larXiv:11 07.0585 
■*! E.B. Sonin, arXiv: 1107.3378 
''^ H. Zhang, C. Lazo, S. Bliigel, S. Heinze, and Y. Mokrousov, 

Phys. Rev. Lett. 108, 056802 (2012). 
^^ S. Heinze, Appl. Phys. A 85, 407 (2006). 



